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Abstract: A simple class of chaotic systems in a random 
environment is considered and the fluctuation theorem is 
extended under the assumption of reversibility. 



1. Random chaos 

A chaotic system in random environment is defined by 

(1) an environment which can be in states ui belonging to 
a space CI and varying in time; the time evolution being 
given by a map t on CI; 

(2) a r-invariant probability distribution Ponfi describ- 
ing the statistical properties of the environment evolu- 
tion; 

(3) a family of maps x — ► fui(x) of a manifold T (phase 
space) into itself. 

The space of the random events Cl will be supposed 
a space of sequences of finitely many symbols u> = 
{ujj}°°^_ 00 , for instance a sequence of spins ujj = ±1, 

and 



with the usual metric = 
with t the shift to the left: (tcj)j = 



+i. A "reflec- 
tion" operator will be defined on Cl as u>< — >uj t , with 



(« 3 



OJ- 



The probability distribution P will be 



a t— invariant and reflection-invariant mixing process on 
Cl, for instance a Bernoulli shift or a Markov process or 
a general Gibbs distribution on Cl. Reflection-invariant 
means that P(E) = P(E T ) for every Borel set E, where 
E T is the image of E under the reflection ux — >u) T . The 
manifold T will be a torus T m and f u will be supposed 
to be close to a map / which is a "linear hyperbolic tor- 
sion" of T , independent of the environment, i.e. close 
to a map which is defined by the action of a matrix with 
integer entries and determinant 1 on the torus. 

For instance the point x € T can be visualized as a 
pair of clock arms x = tp = (<-px,f2) frantically moving 
as 



if -> S (ip) = {(pi + f2,<Pl)- 



(1.1) 



The map is a "linear hyperbolic torsion" of T and is a 
paradigmatic example of Anosov maps. 

The restriction on the choice of /, however, is not dic- 
tated by a simplicity requirement alone. The present 
technique is very special for the cases that we shall con- 
sider. However it should be possible to extend them to a 
more general situation using the ideas in [1]. 

The actual f w will be small perturbations of /. For 



simplicity, for a small e, we will take it of the form 

x^ f u (x) = f(x)+e^(x), (1-2) 

where luq € { — 1,1} and tp is a trigonometric polynomial 
periodic on T m . Therefore we are led to consider the 
dynamical system (Cl x T, S e ), with 



S e (u},x) = (tu, fw(x)). 



(1.3) 



In other words at every instant t the next coin is flipped 
(or the next spin state is observed) and the system point 
xt is moved by f(x t ) + euj t ifj(xt). We shall denote x = 
(u>, x) the points in Cl x T . 

According to standard terminology if /l*o denotes the 
normalized volume measure on T ("Liouvillc measure") 
the system will be said to possess a well defined statistics 
if, for almost all (uj, x) E fix J 7 in the P x /i -distribution 
sense, the limit 



T-l 



lim 



1 

3=0 



F(tJ(*)) 



Kdy)F(y) 



(1.4) 



exists for all continuous "observables" F on T (hence for 
many others), and is independent of x = (lj,x), thus 
defining a probability distribution /i which will be called 
the statistics of the motion or the SRB distribution; here 
ft? ( x ) j j € Z is defined as the composition of j maps 
Us-**, ° ■■■ % if j > and o . . . o f-\ u if j < 0, 
and = 1. More ambitiously one could look for an 
5^-invariant distribution of the form 



^ srb (do;dx) = P(du))^(dx), 



(1.5) 



with limT~,oo ^Y^J=o F(S^(u>,x)) existing for all con 
observables F on 



tinuous observables F on CI x !F, apart from a Im- 
probability set with respect to P x fj, , and having the 

form IaxJ^ F i ujl ' ^Wb(dw'dy) d = (-F) 8rb ; when existing, 
/tx sr b is also called the SRB distribution and it is related 
to (1.4) by n(dx) = J n fi M (dx)P(du). 

Systems more general than the above have been con- 
sidered in [2] where 



cr u (x) = - log 



det 



Ox 



(1.6) 



is introduced and called entropy production rate: further- 
more it is proved that a + = f (cr) srb > 0, see (1.5). Here 
we shall mainly consider the case cr+ > 0, which is the 
generic case. 

Particular attention will be given to reversible systems: 
namely systems for which there is a smooth map X : 

T(u),x) — (<*>' (a;), x'(u>, x)), with uj'(u>) = f uj t , such that 
I 2 = ±1 and 



(1.7) 



for some integer k. 



3. Decoupling and shadowing 



2 



Examples of such systems are rather simple to give in 
the non-random case: however in the random cases (like 
the ones considered here) the very existence of X is not 
obvious. 

In this paper we follow the general ideas of the "cluster 
expansion" techniques which were introduced in [3-7] but 
developed following the tree expansion method of [1] , see 
also [8]. 

2. Fluctuation theorem 

In the non-random case, i.e. if f M = f is independent of 
u>, and for / a general Anosov map it has been shown 
that a time reversal symmetry can be translated into 
certain relations between the probabilities of the "large 
fluctuations" of the time averages of the dimcnsionlcss 
observable a w (x) /a+ . Namely consider the observable 

1 (T-\- 

and call ttt(p G A) the probability that it takes value 
in an interval A evaluated in the stationary SRB distri- 
bution (which for Anosov maps exists). Then itt(p G 
A) = e T C(p)+°( 1 ) where is an analytic function of p 
defined in the interval (—p*,p*), for some p* > 1. The 
time reversal symmetry implies 

C(-P) = C(P)-P*+, \P\<P*, (2.2) 

which is called the fluctuation theorem, and was proven 
in [9-11]. 

The above statements concern non-random maps (i.e. 
f u = f is u> independent). This note is dedicated to 
show that (2.2) also holds for the random noise cases in- 
troduced in (1.3) with f u given by (1.2) and / = So, 
as given for instance by (1.1), when a time reversal 
symmetry, (1.7), I(ui,x) = (u>'(u>), x'(u>, x) holds with 

u)'(oj) = w T = {uJ-j}'jL_ 00 and P is invariant for the 
map u>< — >u; T . 

A proof can be achieved by trying to fall back on the 
already existing proofs of extensions of the fluctuation 
theorem to stochastic processes; the first one in [12] has 
been followed and widely extended in [13, 14]. 

In the first examples mentioned at the beginning of 
Sec.l, in which the u>'s are a sequence of ±1 with inde- 
pendent distribution or with finite range coupling, it is 
clear that the process ^(dwdx) is a Markov process and 
the analysis in Sec. 2 of [13] can in principle be applied. In 
the present case the phase space is a continuum but the 
results in the quoted paper nevertheless apply, at least if 
one is willing to consider formally ratios of delta functions 
and interpret them as suitable Jacobian determinants. 

To discuss the example of a general Gibbs distribution, 
in which the distribution of the u> is far from Markovian, 
along the lines of [13] is harder. But one can consider the 
more general approach in [14]. However in this case one 



has to prove that the SRB distribution is a "space-time 
Gibbs state", which is the basic object studied in [14]. 

In all cases an explicit determination of the station- 
ary state seems necessary in order to check the assump- 
tions and to define and compute the quantity e(A) of the 
quoted references. In other words the existing literature 
provides strong arguments (particularly in the Marko- 
vian cases) for the validity of a relation like (2.1): but 
some work to check it remains to be done, on a case by 
case basis, and a substantial further work is necessary if 
one wishes to compute the function ((p) or at least some 
of its properties. 

We restrict attention to the (very special) cases 

(1) S £ defined by (1.2), (1.3), (1.1) and 

(2) ip( x ) trigonometric polynomial. 

We shall show that the technique of [15] can be applied 
to obtain, for e small enough, the results: 

(1) the SRB distribution \x exists and is unique; 

(ii) it is a space-time Gibbs distribution; 

(iii) cr+ > at least for £ ^ small enough; 

(iv) under the extra assumption of existence of a time 
reversal symmetry for S, the large deviation rate £(p) for 
p(x), see (2.1), satisfies the symmetry relation (2.2), i.e. 
a fluctuation theorem holds. 

Remarks: (i) The extension of the results to perturba- 
tions of torsions of tori of arbitrary dimension is not 
treated because it is identical to the one we present here 
but the notational burden hides the simplicity of the 
ideas. 

(2) The apparent restriction to perturbations of Anosov 
maps generated by a linear torsion of a torus is due to 
the use made of the flatness and parallelism properties 
of the stable (respectively unstable) planes for such tor- 
sions. However by applying the methods in [1] the results 
can be extended to the general case in which perturba- 
tions of generic analytic Anosov maps of tori of arbitrary 
dimension are considered, see comments in Sec. 8. 

3. Decoupling and shadowing 

The key for fulfilling the program set in Sec. 2 is to show 
that there is a map x' = (u)',x')< — >x = (u>,x), that we 
call H £ , of the form 

{u',a/) = H e {u,x)^(u,x + h u {x)), (3.1) 

with h M (x) analytic in e (but not in x, in general) such 
that 

S e oH e = H £ oS Q , (3.2) 

where Sq(u),x) = (tuj,Sqx). This means that there is 
a change of variables turning the perturbed map into 
the unperturbed one, i.e. the perturbed map S £ can be 
"conjugated" to the unperturbed one So- The relation 
(3.2) becomes, see (1.2), 

Soh w (x) — h TM (Sox) = —eip w (x + h U} (x)), (3.3) 
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where So is a 2 x 2 matrix in our case. 

We look for a solution which is analytic in e: h u! (x) = 

ehu\x) + e 2 hu\x) + with h^\x) e-independent 
functions. For instance the equation for the first order is 

S h£\x) - h<U(Sox) = -1> u {x). (3.4) 

Call v± the two normalized eigenvectors of S relative to 
the eigenvalues X± and let A the inverse of the largest 
one (A = |(V^-1)), so that A+ = A _1 ,A_ = -A: A < 1. 

The functions ip u , h M can be split into two components 
along the vectors v± : 



tp w (x) = ip Wi+ (x)v + + il> Wi -(x)v-, 



(3.5) 



h M (x) = h Ml +(x)v+ + h Wl -(x)v-, 
and the equation (3.5) for h$± gives 

\-h$_(x) - h^,_(Sox) = ^,-(x), 
which can be solved by simply setting 



(3.6) 



where a = ±, Z+ = [0,oo)nZ, Z_ = (-oo,0)nZ and 
the inequality A < 1 ensures convergence. 

(k) 

Hence the equations for h± become 

ew-Ei e • 

s = feiH |-fcs = fc-l, fcj>0 

"1 «.=± 

s 

• E « A « b+1| "(llK- • 0*))VVp«,a(Sfc)- (3-8) 

■ (n^ >a ,(«). 

i=i 

This can be written via a graphical representation as in 
Fig.l, where the "graph element" on the l.h.s. represents 
hu]a (x) ■ Representing again, in the same way, the graph 
elements that appear on the r.h.s. one obtains an expres- 
sion for h^ a (x) in terms of trees, oriented "toward the 
root". 

(hi) 
(to) 




(k s -i) 



Fig.l Graphical interpretation of (3.8) for k > 1. 



A tree 9 with k nodes will carry on the branches I a 
pair of labels ag,p£, with p<eZ and ai G {— , +}, and on 
the nodes v a pair of labels (not shown in Fig.l) a v ,p v , 
with a v = ae v and p v £ Z Qij such that 



P(v) = Pi v = ^Pv, 



(3.9) 



where the sum is over the nodes following v (i.e. over 
the nodes along the path connecting v to the root), £ v 
denotes the branch exiting from the node v, and to each 
tree we shall assign, given w, a value given by 



Val(*)= J] -^A^+ 1 I- 

Sy ! 



vev($) 



(f[d aVj )^ MM , av (S p Q ^x), 



(3.10) 



where d a • d x , V{9) is the set of nodes in 9, the 

nodes vi,...,v 8v are the s„ nodes preceding v (if u is 
a top node then the derivatives are simply missing). If 
Qk,a denotes the set of all trees with k nodes and with 
label a associated with the root line, then one has 



) = E £fc E Vai w> ( 3 - u ) 



fe=i 



oee k , 



and the "only" problem left is to estimate the radius of 
convergence of the above formal power series. For this 
purpose it is convenient to study the Fourier transform 
of the function h Uja (x). This is easily done graphically 
because it is enough to attach a label v v e Z 2 to each 
node and define the momentum that "flows" on the tree 
branch £ v as v lv = J2 w ^v u w 
Then (3.11) becomes 



(3.12) 



k=l 



with 



h {k) = V 
0ee fe ,„, o 



vev(e) 



(3.13) 



where Qk,v.a denotes the set of all trees with k nodes and 
with labels v and a associated with the root line, and Vo 
is the node immediatley preceding the root. 

Since we are assuming that ip u is a trigonometric poly- 
nomial of degree TV and with uniformly bounded coeffi- 
cients as u> varies we can define F = max„ IV'w^l- Then, 
to study the regularity of the conjugation, we shall esti- 
mate ^2 V \v\P\h£),\ for P > 0. The only problem is given 
by the presence of the factor \v\& . In fact consider first 
the case (3 = 0: there are only (27V+1) 2 < (3iV) 2 possible 

choices for each u v , given p(v), such that IS^^Vy] < N. 
Hence fixed 9, {a v } v( zv(0) an d {Pv}vev(0) the remaining 
sum of products in (3.13) is bounded by (if A = A+ 1 = 
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-A-) 



(3N) 2k N k F k ] [ 



vev(0) 



(3.14) 



having bounded by N k the last product in (3.13). The 
sum over the {p v } v ev(0) i s a multiple geometric series 
bounded by (2/(1 - A)) fe . 

The combinatorial factors arising from the Taylor ex- 
pansion leading to (3.8) (see the s! _1 factors in (3.8) or 
the s^r 1 factors in (3.14)) compensate the large num- 
ber of trees so that the factor Ilt;(l/ S « becomes, after 
summing over all the trees, simply bounded by 2 3k , (2 2k 
due to a (poor) estimate of the number of trees, after 
appropriate combinatorial considerations, and 2 k due to 
the sum over the labels a v = ±). 

Therefore for (3 = the conjugating function H £ exists 
and it is uniformly continuous and uniformly bounded 
with a uniformly summable Fourier transform, for e (even 

complex) in |e| < e (0) d = (3jV) _3 -F _1 2~ 4 (l - A), where 
an extra factor 2 _1 has been inserted in order to obtain 
uniform bounds. 

Taking (3 > requires estimating we bound it 

by E„ Thcn from l-V^^I ^ N (otherwise the 

value vanishes), we infer that \u v \ < X'^^BN, where 
B > 1 is a suitable constant. The sum ^2 V \v v \P is over k 
terms which can be estimated separately so that we can 
write X)ul iy «l /3 — where \vy\ = m.ax. v \v v \. This 

can be taken into account by multiplying (3.14) by an 

extra factor (BN)^X~^ P ^ < BNX~ ^« M . There- 
fore if (3 < 1 the Y, v W\ fi \ h *U is bounded as in the case 
(3 = but the result is modified into 



e (/?) = (SN^F-^l - A 1 - /3 )2- 



(3.15) 



This shows that H e is analytic in e in the disk with ra- 
dius Eo{f3). Furthermore, since in (3.15) we inserted (for 
simplicity) an extra factor 2 _1 in excess of the result ob- 
tained by the procedure described, the Holder modulus 
is also uniformly bounded by a suitable function C{(3) of 
(3. Note that e {(3) -> for (3 -» 1. 

The map H e is a homeomorphism of T 2 at fixed ui. 

Remark: The result can also be interpreted as a "shadow- 
ing theorem" . Consider the "noisy" trajectory k — > S k x 
as a perturbation of a noiseless one. Then a noisy tra- 
jectory starting at x = (u>, x) will remain forever close 
to the trajectory of the point iJ~ 1 (w,a:) evolving under 
the noiseless motion. This is usually called a shadowing 
property. 



(a) at every point x = (lj,x) e il x T the plane W(x) 
tangent to T contains two planes, W^(x),W%(x) with 
W(x) — W£(x)®W£(x), with positive dimensions d s ,d u 
and which are covariant, i.e. SW%(x) — W" ul (f u ,(x)); 

(b) for some C, n > 0: 



\\d x f™(x)v\\ <Ce 
*(-k) 



-kn I 



v € Wl(x), 



d x fT k \x)v\\ < Ce- kK \\v\\, v e WZ{x), 



(4.1) 



for all k > 0; f^ k is defined below (1.4) and the sizes 1 1 • 1 1 
of the vectors are evaluated in the metrics at the points 
to which they are applied (i.e. f^ k (x) or x or fv~ k \x)); 
(c) W®(x) are continuous in u:,x. 

In the example (1.3) with e = one has d x (jt} ±k ^) — 
SQ k so that properties (a,b,c) are satisfied, trivially, be- 
cause the map considered is an Anosov map. 

It is also natural to define the expansion and contrac- 
tion rates at (u:,x). They will simply be 



log 



dct 



dfJ aN \ x) 
dx 



) , a = u,s, (4.2) 

/ a 



where (aJV) means N if a = u and — N if a = s while 
the label a appended to the Jacobian matrix means that 
it is regarded as acting as a map from the plane W®(x) 
to the plane W? N Jf* N (x)). 

Note that the above expansion rates are not intrinsic 
geometric objects as they depend on the metric used on 
the manifold, as also the constants C, k do. 

The volume measure /xo(dx) on T can be visualized, 
fixed w S fi, in terms of the above rates: imagine to 
fix, for every point x G J 7 , two segments located with 
center in x and lying on the unstable or on the stable 
planes through x, and call them A™ (a;) and A^(x), re- 
spectively. Then we can build a tiny rectangle through x 
by considering the line elements through x 



S s (N,x) = f^ NM (A 



'w(/w 



*(-N) 



(*))), 



(4.3) 



and by drawing through each point y £ S s the line ele- 
ment S w (N, x) and through each point of y € S s the line 
element 6 s (N, x) and collecting the intersections: this is 
well defined, and will be called S u>i n(x). This is close to 
a rectangle with volume exponentially small as N — > oo 
and, to leading order in N, given by 



vo!5 u ,n(x) 



; g JZ,N( X ) 



(4.4) 



with, if a w (x) equals the angle between the stable and 
unstable planes at x, 



4. Heuristic considerations 

The notion of hyperbolicity or chaoticity can be ex- 
tended, [16], to random systems, for a general map S 
acting on ft x T, as follows: 



lul {x) A ^ cosa^x) J] \A a T(aN) Jf^ N Hx))l (4.5) 

a—u,s 

where (aJV) has the meaning explained after (4.2). The 
function J w (x) is uniformly bounded away from and 
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oo. Hence the ratio of the volumes of two such paral- 
lelepipeds, centered at x and y, can be computed by the 
ratio between the contraction factors in (4.3) evaluated 
at the two points (because the basic surface elements 
A£(x) have all the same size and shape up to a factor 
1u{x)h u {y)). 

A consequence is that if we attribute to each such par- 
allelepiped a measure proportional to 

e -Ja,,2Jv(/<- *) 5 (4.6) 

then we expect that in the limit as N — > oo the 
probability distribution has a limit ^(dx) such that 

fi M (dx)P(du:) = f fi(dxduj) is an invariant distribution for 
the system (O x S £ ) which should describe the statis- 
tical properties of almost all data initially chosen with 
the distribution no(dx)P(duj): i.e. fj,(dxdu>) should be 
the SRB distribution. Likewise replacing in (4.6) N with 
— N and J u with — J s one should obtain the SRB distri- 
bution for the backwards motion (i.e. for the map S^ 1 ). 

5. Overshadowing 

Therefore we look for an explicit algorithm to construct 
a useful representation of (4.6) or, what is the same, for 
the function J M ,i{x) because 

j=-N 

by the chain differentiation rule for composition of func- 
tions and by the multiplication rules of determinants of 
matrices. 

We shall use the notations of Sec.3, where the con- 
jugation H £ , transforming a perturbation of the ma^ 
S e {w,x) = (t«, So x + eip w (x)) of the torus fl x T 
into a noiseless map So(u>,x) = {tuj,Sox), has been de- 
rived. The homeomorphism (u>',x')< — x + h M (x)) = 
H e (u>, x) can be used not only to construct the dynamics 
but also the stable and unstable manifolds of each point. 
The latter manifolds through (u>, x + h M {x)) are given by 
parametric equations of the form 

1 a {t) = {u,h u (x + tv a )) teR, a = ±, (5.2) 

where t — > x + tv a is the unstable or stable manifold for 
the unperturbed map So of T 2 into itself, respectively if 
a = + or a = — . 

The construction of Sec.3 gives a Holder continuous 
conjugation with a prefixed exponent j3 < 1 for a per- 
turbation strength e that is suitably small, depending on 
how close (3 is to 1. However, in general, the conjugation 
is not be diffcrentiable: therefore (5.2) cannot be used 
to compute the derivatives appearing in the Jacobians 
Jui i( x ) an d the parametrization (5.2) is not very useful. 

Instead of constructing the stable an unstable mani- 
folds parameterized so that the required derivatives ap- 
pearing in (5.1) can be computed (or just shown to exist) 



we remark that all we need are the expansion coefficients 
in the stable and unstable directions of any point (u>, x). 

To find them we apply once more the technique dis- 
cussed in the previous two sections. Calling the {non- 
compact) space OxT xl we define on it the dynamical 
system (6, So) with 

5o(w, x, v) = (tw, S x, S v), (5.3) 

where v has to be interpreted as a tangent vector to T 2 
in x and So is the action of the derivative of the map 
x — ► Sqx on the tangent vector v. Since the map So is 
locally linear there is no point in distinguishing between 
the constant matrix 8 x SqX and the map So so that we 
indulge in the abuse of notation implicit in (5.3). 
We can consider the perturbation of So defined by 



= (ra, S x + eipuj(x), S v + e(v ■ d x )ip w (x)), 

where the third term is the vector v tranformed by the 
map x' = x + Eip u (x) into a tangent vector at x'. 

For uniformity of notation in what follows, we shall 
denote Suj(x) = Sqx + ei/) M (x) instead of f M (x) as done 
so far. 

This system fails to be an Anosov system not only be- 
cause of the noise (which acts trivially, however) but also 
because the space fi x T 2 x R 2 is not compact. Never- 
theless we can still try to find an isomorphism between 
S E and Sq- 

Success would mean that it is possible to decouple from 
the noise also the evolution of the tangent vectors, i.e. of 
the infinitesimal displacements: not only we could match 
individual trajectories in the perturbed and in the un- 
perturbed system, as done in Sec.3, but we could even 
achieve "shadowing" of infinitesimally close pairs of tra- 
jectories which would split apart at the same rates as 
the unperturbed ones. If possible this property could be 
called overshadowing: but not surprisingly this turns out 
to be in general impossible (as it will be implicit in what 
follows). The next simplest goal, in trying to compare 
infinitesimally close pairs of trajectories, is to try to con- 
jugate S e with a map Sq, e defined by 

<So, £ (a>, x, v) = (tw, S a x, S a v + T u) (x)v), (5.5) 

with T u} (x) a 2 x 2-matrix commuting with So (hence 
diagonal on the basis v±), analytic in e and vanishing at 
e = 0. Therefore we look for a map H e of a simple form 
and such that S e o H e = H £ o So, e , i.e. 

H £ : (u;,x,w)< — ► (u>, x + h M (x)), w + K w (x)w), (5.6) 

where K is a 2 x 2-matrix analytic in e and at e = 0, 
to be determined. In guessing (5.6), advantage is taken 
of the already known conjugation H e between S E and So, 
from the analysis of Sec.3, giving us the function h ul {x). 
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Remarks: (1) Let fC^(x) = 1 + K LJ (x) and let C ul {x) — 
So + T u> (x) be functions with values in the 2x2 matrices. 
Then the conjugation (5.6) is equivalent to the following 
equation 

d x S u {H e {(jj,x))K u {x)v = JC TU ,(S x)£ w (x)v (5.7) 

and to (3.3) for H e , which is an equation independent 
of (5.7) (already solved in Sect. 3). This implies that the 
vector w Mt ±(x) = 1C u (x)v± satisfies: 



d x S u {H u {x))w Ut ±{x) = X U! ±(x)w TU1 ,±(S x) = 
= (A± +'f±,u>{x))w TUJ ,±{S x), 



(5.8) 



where X^.±(x) are the diagonal elements of £^, £ (x). 

(2) The conjugation defined in (5.6) is the right one to 
look at the stable and unstable directions. Its solution 
gives the stable and unstable directions as functions of 
the unperturbed point i/ 1 T 1 (a;,a;). 

(3) Equation (5.7) does not determine K^x) uniquely. 
Indeed, if lu>,±(x) are two non-zero functions from T to 
K and if \u>,±(x), w Mt ±(x) solve (5.8), also 



w w ,±(x) = L,±(x)w M ,±(x) 



(5.9) 



solve it. To fix this ambiguity we will require that the 
diagonal elements of K lA) (x), on the basis v±, are equal 
to 0, i.e. the matrix K w (x) is completely off-diagonal. 

To simplify the rather heavy notations we denote 

ip a {u),x) = V>u>,a(z), H(u),x) = (u),x + h w (x)), 

T(u:,x) = T M (x), K(u>,x) = K M (x), 

e = x=(w,a;), (5.10) 

and the equation that the matrices have to 

satisfy is a transcription of (5.7) 

= -ed Xj iPi(H(0) - ed x ,il>i(H(£))K(Z) aj + (5.11) 

where d x denotes a derivative of ip with respect to its orig- 
inal argument and repeated indices mean implicit sum- 
mation (to abridge notations). 

We write the above matrix equation on the basis in 
which So and T are diagonal, i.e. on the basis formed 
by the two eigenvectors v± of So in which the matrices 
K u , r u have been assumed to take the form 



r(0 = 



7+(0 






7-(0 



o MO 

MO 



(5.12) 



If a = ±, (3 



-a and a n — v n 







d v , (5.11) becomes for 
ed a MH(0) - eKf, ia (Qdf,il> a (H{0)+ 



+ 7a(0. 



and, for i ^ j and if A + = A , A_ 
values of S (with A = (y/E- l)/2), 



(5.13) 
-A are the eigen- 



- ed^ a (H(0) - sK a ^)d a i> a (H(x))+ (5.14) 

+ K aJ ,(So(0)MO- 

If a' = SL ^, and a" — we can rewrite the equa- 
tions (5.13) and (5.14) as 

j a (O=ed a ^ a (H(O)+ek (Od ^ a (H^)), (5.15) 

and, respectively, 

MO + ^M^O = aX(- ed ^ a (H(S^))- (5.16) 

- ek a (SZ'Od a MH(SZ'0) + k a (sfC)lfj(SoO)- 

These equations are in a form suitable for a recursive 
solution in powers of e. For instance the first order is 

fc^CO + X 2 k^(SoO = -Ad_V+(0, (5-17) 
fc a) (0 + aV^^O = Xd + ^(S -H), 
which has the solution 

l { a\0 = d a M0 « = ± 

oo 

fc + ) (0 = -AE(- 1 ) 9A295 -^+( 5 oO, 



(5.18) 



9 =0 

OO 



fc a) (0 - a£(-i)«a 2 ^ + v>-GS - (<z+1) 0- 

9=0 

The equations for 7^(0 and fc Q fc ^(£) can be represented 
in graph form by suitably modifying the similar represen- 
tation derived for h!$ a {x) in Sec. 3; compare Fig. 2 with 
Fig.l. 

(fe) (fe-i) >a 



(fe) 



/3a 

oVW» - 

(fe-1) 



^(fe-l-p) 

(P) _a_J^ P {P) 



Fig. 2 Here a = ± and /3 = —a. All the lines have to imagined 
to carry arrows (not drawn) pointing toward the root. The line 
carrying a label a and emerging from a circle or a square with label 
(fe) denotes 7^ or , respectively. The wavy line emerging from 
a bullet with label (p), with 1 < p < fe — 1, ending in a small circle 
and carrying a pair of labels 7,5, represents [dyips]^ (H(£)), the 
p-th order in the power expansion in e of (d^^i) (if(£)). The 
small square in the node closest to the root into which a wavy line 
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arrives will carry a label q = 0, 1, 2 . . . (see the first two graphs 
in the second line) and it expresses that H(S£ ( in thc 

second line), with q+ = q and q_ = q+1, is thc argument in which 
the functions djipg are computed. Furthermore a summation over 
q = 0, 1, . . . and a multiplication by — q( — l) <3 A 1 + 2<? is understood 
to be performed over the nodes represented as small squares. The 
last graph in thc second line docs not carry any further labels and 
it represents the last term in (5.16). 

The representation is drawn in figure Fig. 2 and the 
symbols are explained in the corresponding caption: the 
reader will recognize in them a pictorial rewriting of 
(5.16). 

In this case too we can continue the expansion until in 
thc r.h.s. of Fig. 2 all top nodes of the graph are either 
squares or circles carrying a label (1), i.e. they represent 
a first order contribution to T (circle) or to K (square), or 

bullets representing cither 9 Q V/3(£) or d a ^p/3(Sg +Q9a £); 
cf. thc caption of Fig. 2 for notations. Of course the 
latter quantities can themselves be represented by the 
tree expansion discussed in Sec. 3: if we do so then we 
obtain a full expansion in powers of s in which the wavy 
lines with label p are replaced by a tree with p nodes. 

The rule to construct the value of each tree graph 
is easily read from (5.11) and from thc rules discussed 
above, and in Sec. 3, to build thc value of trees represent- 
ing h. 

The estimate of the fc-th order contribution is given by 
(3.14) with an extra factor TV* to take into account the 
extra derivatives due to the lines with two labels. Also 
the counting of the trees has to be modified but the end 
result will be that T and K are expressed by convergent 
series in e for |e| < So(f3), where So(/3) can be taken of the 
form (3.15) with a different numerical factor and with N 
replaced by N 2 . 

6. Construction of the SRB distribution 

Consider the partition £ of O x T by the sets E^^ = C u x 
P a where C — {C\, . . . , C m } is the partition of O accord- 
ing to the value of the symbol uo and Vq = {P\, . . . , P n } 
is a Markovian pavement of T for the map So, [8]. It is 
convenient to suppose that the Markovian pavement Vo 
is "reversible". This means that there exists a map Io, 
with the property that Iq = ±1 and Io°Sq = SQ k °Io 
for some integer k, such that for all a — 1, . . . , n one has 
IoP a = P a T for some a T (this is not restrictive as it can 
be achieved by considering the new pavement loPo^'Po)- 
Given a point £ let (<*>, cr) be its "history" on the parti- 
tion £ for the map <Srj. This means that Sq£ e E LUky(Jk for 
all k. In terms of the sequences (u>, cr) the motion is sim- 
ply a translation (r,p) such that (ru:,pcr)j = (u>,er)j + i, 
[8]. 

The partition £' = H e {£) is then Markovian for the 
perturbed system. A correspondence can be generated 
between points in O x T by associating two points £' and 
£ if the first has history (u>, cr) on £' identical to the his- 
tory of £ on the partition £ under the action of the map 



So- we denote £' = X e (u;,cr) and £ = Xo(uj,ct). The 
correspondence so defined is identical to the correspon- 
dence between £' and £ defined by £' = H(£). And in 
Sec. 5 it has been proved that 

There exists £q{/3) such that the expansion rate A Uj0l ,(<t) 
of S £ along the unstable manifold of £' is defined and 
holomorphic in e in the disk \e\ < £o(P)- As a function 
of (w, cr) it is Holder continuous with some exponent < 
(3 < 1 and modulus C({3). 

This can be seen as follows: one just notes that if 
£' = H e (£) then thc unstable direction at which is 
a vector tangent to T ', will be w UJi+ (a;) = v + + K w (x)v + . 
Then the expansion rate along the unstable manifold of 

C = {u,a>) = H e (£) will be e A »^ def 
by (5.8), (4.1), (4.2), 



|J2J,i(^)| and, 



A^.uicr) = log ((A+ + j u ,+(x)) 



\Wu,+ {S x)\ 

\w w ,+ {x)\ 



), (6-1) 



where |w Wi+ (a;)| = y/1 + k W}+ {x) 2 . 

The functions 7o>.+ ,&a>.+ are analytic in e and Holder 
continuous in £ with a uniformly bounded modulus C((3) 
if (3 < 1 and e is small enough. Since {(uii, a i )}°l_ oc 
fixed means £ fixed (because cr is the history of a point x 
on a Markov pavement for the e-independent map So), 
we see that A^ tU {cr) is analytic in e at fixed (w,cr) and 
Holder continuous in cr) and therefore in £. 

Hence the function A u , tU (cr) generates a short range 
potential on the one-dimensional "spin system" whose 
states are the sequences (u), cr), in the sense of [14]. 

An identical argument holds for the similar function 

A u} . s {cr) A = log \Ju\{x?) | ■ It follows, by standard ar- 
guments on onc-dimensional Gibbs distributions, see 
Chap. 6 in [8], that the "short range" Gibbs distribution 
v with formal energy function 



E 



-4 T 



k=- 



i (p k cr)+Y / A T 

k=0 



,u{p k(T ) 



(6.2) 



is well defined and it has the property that if C^ CT is 
defined as {a/, cr' \ uo[ = Ui, u\ = Ui, V|i| < N}, then 



B- 1 < 



vol(C£ CT ) 



< B 



(6.3) 



for a suitable constant B > 0. 

Therefore the distribution v is absolutely continuous 
with respect to the volume distribution (since j w in (4.6) 
is bounded above and below). 

The distribution v is not translation-invariant because 
its Gibbs potential in the "far future" is governed by 
A^. u {cr) while in the far past it is governed by j4 WjS (<t), 
see (6.2). But short range Gibbs distributions (whether 
invariant or not) enjoy exponentially mixing properties: 
hence the distribution v, and consequently the volume 
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distribution, will be very close on the evolution of the 
cylinders, for very large times, to the invariant Gibbs 
distributions with (formal) potentials 

E. , k \ as v — > +00 it a = u, ,„ AS 
A ^AP k ^ {asJ--ooifa = a ; (6 - 4) 

k— — 00 

respectively. Therefore the distribution with potential in 
(6.4) with a = u is the SRB distribution for the forward 
evolution while for a = s the distribution is the SRB 
distribution for the backward evolution. 

7. Time reversal symmetry 

The result of Sec. 6 allows us to conclude that the SRB 
distribution /i sr b will be a "space-time" Gibbs state for 
the energy function A UtU {(r) in the sense of [14]. 

Time reversal symmetry is a symmetry that should be 
inherited from the microscopic time reversal symmetry 
of purely Hamiltonian dynamics. It may be destroyed 
by the thermostats, or better, by the models that are 
introduced to describe them. 

In presence of noise the matter is even more involved. 
However under the assumptions that we have posed upon 
S e it is possible to show that if the noiseless system So is 
reversible a "kind" of time reversal symmetry will con- 
tinue to hold when the noise is not neglected provided 
the noise also has a reversibility property. 

We supposed, since the very beginning in Sec.l, re- 
versibility of the stationary probability distribution P on 
in the sense that the operation 9 defined by 9lj = u> T 
with wj = oj-i ("history reversal") conserves the distri- 
bution P: 

P(E T ) = P(E), (7.1) 

where E T = {ui T : u> <= E}, for all (Borel-)sets E. 

Suppose that the unperturbed system evolution on T ', 
that we assume to be an Anosov map /, is reversible 
in the sense that there is a map Iq : T — > T with the 
property Io°.f k = f~ k °Io f° r some k and J5 2 = ±1, or 
more generally Iq =1 for some k' even. 

The above assumption is satisfied in the examples 
treated in the previous sections because the time rever- 
sal operation I ((pi 7 (p2) = (V2>— Vi) nas the property 
I So = S _2 °/ , if S is given by (1.1). 

The problem is that it is not clear that a time reversal 
symmetry holds when there is a perturbation, even if 
small. But, if e is small, the time reversal symmetry 
Zo(a>, x) = (co T , Iqx) is not "broken" by the perturbation 
and the noisy system admits a time reversal symmetry 
I £ . The latter can be immediately written as 

l E =H s oI oH-\ (7.2) 

However for the same reasons for which the conjugacy 
between S e and Sq was not smoother than Holder con- 
tinuous it will turn out that X e is also not necessarily 
smoother than Holder continuous. 



In terms of symbolic dynamic the time reversal X e acts 
on a point whose symbolic history is {us, a) by trans- 
forming it into the point with symbolic history (u> T , cr T ), 
where (f T )i = c 7 ^, if time reversal transforms the ele- 
ment P a of the Markov pavement for So into P a T . It is 
therefore clear that = 1 if Iq =1. 

The lack of smoothness of the new time reversal im- 
plies that the fluctuation theorem can certainly be proved 
only in cases in which the perturbed system also admits 
a smooth time reversal. Then the proof is the same as 
the one for Anosov maps because the latter was based 
on the existence of a Markovian partition of phase space 
and on the Gibbs property of the SRB distribution with a 
Gibbs weight generated by a function A UtU} (cr), as shown 
in Sec. 6, satisfying A umj t(ct t ) = —A s ^(cr) (which fol- 
lows from time reversal symmetry and from the isometric 
nature of such reversal) [9, 11]. 

Examples can be obtained easily: consider the paradig- 
matic example S e of the previous sections and the map 
T(w, x, t], y) defined on Q x T 2 x Q x T 2 by 

T(u,x,r,,y) - (S e (u, x), ^(77, y)). (7.3) 

Then a time reversal symmetry (i.e. an isometry anti- 
commuting with time and squaring to the identity) can 
be defined as X((u>, 77, x, y) = (rj, uj, y, x) and the fluctu- 
ation theorem applies to this system. 

The last example is somewhat artificial: in applications 
time reversal should be a built-in symmetry so that its 
checking should be immediate. This is often the case in 
non-random systems. Note that in numerical simulations 
realizing the system T out of a system S would be easy. 

8. Structural stability 

All what has been said about the existence of the conju- 
gacies to the unperturbed system can be extended to the 
cases in which the torus has arbitrary dimension and So 
is an hyperbolic linear torsion (i.e. it is a map generated 
by a matrix So with determinant ±1, integer entries and 
all eigenvalues with modulus different from 1), by replac- 
ing A± by the determinants of So restricted to the stable 
or unstable planes. 

Consider the case in which the map x — > f(x) is just 
an Anosov map of a general manifold J- ', still assuming 
analytic regularity. If the unperturbed Anosov map / is 
close to some hyperbolic linear torsion So of a torus it will 
be possible to find an analytic conjugacy to the systems 
to So simply because we can think that the perturbation 
eipu, contains a non-random part which is the perturba- 
tion of So into /. The same can be said of maps that can 
be analytically conjugated to hyperbolic torsions. 

To extend the results on the fluctuation relation, how- 
ever, the strong assumption about time reversal has to 
be always added. 

More interesting and subtle is the further question of 
how much the above analysis depends on the assumption 
that we are perturbing a system close to one So in which 
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a reversible noise and a linear torsion of a torus So do 
not interact. 

Under the only assumption that So is analytic and 
Anosov the equations for H will have to be written in 
coordinates and possibly several coordinate charts will 
have to be used. An apparently unsurmountable diffi- 
culty may seem that the directions v± will now depend 
on x and therefore the solution algorithm will simply fail. 
A more attentive study of the equations shows, never- 



theless, that the equations are soluble order by order as 
at every order the curvature of the manifolds systemati- 
cally gives contributions of higher order, see [1] for details 
showing that the existence of a convergent power series 
in e for H exists. Since that is the main difficulty there is 
little doubt that the analysis that we have performed in 
detail for the case of the random perturbations of (1.1) 
carries unchanged to the last case. 
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